Dynamic response of open cell dry foams 
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We study the mechanical response of an open cell dry foam subjected to periodic forcing using 
experiments and theory. Using the measurements of the static and dynamic stress-strain relation- 
ship, we derive an over-damped model of the foam, as a set of infinitesimal non-linear springs, where 
the damping term depends on the local foam strain. We then analyse the properties of the foam 
when subjected to large amplitudes periodic stresses and determine the conditions for which the 
foam becomes optimally absorbing. 

PACS numbers: 



Foams of different types form an important part of 
our daily lives having diverse applications depending on 
their material properties, e.g. density and flexibility [D- 
01 ■ Thus while open cell polymeric low density foams are 
used for car cushions, semi-rigid foams form packaging 
materials while rigid foams can be engineered to provide 
thermal insulation. An important physical property of 
foams is the efficiency with which energy is dissipated in 
these systems that make them ideal as shock absorbers [H, 
4j. Since material properties are intrinsically connected 
to the underlying microstructure 5[ , current research has 
focussed on its classification @,|7| and an exploration of 
mechanical properties using experiments numerical 
simulations [5|, [l(| and theory 



X-ray microtomography |9( and scanning electron mi- 
croscopy of foams reveal their microstructure as con- 
nected random polyhedra with their bulk mechanical 
properties governed by Gibbs elasticity of the faces for 
closed cells, or nonlinear elasticity of struts for open cell 
foams respectively 11 1. Recently Kraynik et al. 0, 0] 



have performed detailed numerical calculations to classify 
foam microstructure for monodisperse and polydisperse 
foams. These simulations indicate that cells having an 
average number of re 14 faces and re 5 edges are the most 
probable polyhedra occurring in monodisperse foams in 
agreement with experiments by Matzke [14| . The sur- 
face energy of such random foams have been estimated 
using scalingarguments and arc in good agreement with 
simulations [7[. 

Analysis of strut deformation incorporating bending, 
stretching 11, 15|, and torsion [l]| have been employed 
to investi gate the linear and nonlinear elastic behavior 
of foams [12[. The nonlinearity arises as a result of 
the coupling between geometry and strut elasticity. The 
micro-mechanics of foams where constitutive models for 
struts are inherently nonlinear have garnered recent in- 
terest [III . The nonlinear mechanical response has many 
fascinating features including strain localisation 17|, bit , 
crus hing H [III], energy dissipation, and shock mitiga- 
tion [20(. 

A complementary approach to the finite element sim- 
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FIG. 1: Static compression data for foams, stress Pa vs % 
strain. Two polynomial forms (see text) are fitted to the lin- 
ear and non-linear part of the constitutive curve. Top curve, 
dense foam; bottom curve light foam. 



ulations as a mean to connect microscopic structure to 
mechanical pro per ties is modelling foams as a hyperelas- 
tic continuum 



2l[ |22|. Motivated by these studies, 



this Letter we propose and validate a phenomenological 
constitutive equation to model the mechanical response 
of open cell solid foams subject to dynamic loading p3| . 
While preparing this manuscript we came across recent 
work in this area by Del Piero et al. [24||. Our model 
is qualitatively different from prior approaches in having 
a nonlinear constitutive equation and a dynamic friction 
coefficient 7(e) that is dependent upon the local stress of 
the system e. Using this model we find excellent agree- 
ment between experimental data for dynamic loading ex- 
periments for two model foams, for low amplitude, low 
frequency oscillations. This leads to several interesting 
results, namely the identification of a scale factor p that 
can be used to predict system size dependence of foam 
properties within the validity of the continuum approxi- 
mation. We then explore, high amplitude, high frequency 
forcing of model foams theoretically and obtain optimal 
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values of p for energy dissipation and phase dccohercncc. 

We carried out static and dynamic uniaxial compres- 
sion tests for two foams with densities p^ ss 41kg /m 3 and 
pi w 19kg /m 3 (called dense and light respectively). The 
densities were measured using gravimetric method with 
cuboid volumes measured using digital callipers and its 
mass using a precision balance (Precisa aB 320M) with 
an accuracy of 1 milligrams. Static compression tests 
were carried out using an Instron model Extra between 
two steel plates having an upper plate diameter 57.3 mm, 
and a lower plate diameter 150 mm. The tests were car- 
ried out at 50% relative humidity and at a temperature of 
20°C. Cylindrical samples of light and dense foams hav- 
ing diameter d = 37 mm and d = 18.8 mm and heights 
L = 21 mm, and L = 12.4 mm respectively were used. 
A compressive displacement rate of 0.167 mm/s which 
translates to a compressive strain rate of 7 x 10 _3 s _1 . 
The voltage output was fed into a Pico ADC-20 data 
logger and compressive stress-strain curves calculated. 
Dynamic measurements were carried out in compression 
using a Rhcomctrics RSA II solids analyzer. The samples 
were cut to a diameter of 19 mm and a length of 12 mm 
and tested between two steel plates whose diameter was 
25 mm. Initially a static compressive strain between 0.01 
and 0.7 was applied. The system was equilibrated for 1 
min following which a dynamic test with a strain ampli- 
tude of 0.001 and frequency / = 1 Hz was applied. The 
storage and loss modulii E' and E" as well as their rela- 
tive phases tan [8] were calculated for each value of base 
compressive strain (increased in increments of 0.02). 

Fig.[T]shows constitutive relation between the compres- 
sive stress a and strain e for the static light (circles) and 
dense (diamonds) respectively. The experimental data is 
fitted (solid line) to a nonlinear function of the form 



e < e : k e 

e > e : £o + h (e - e ) + K (e - e ) £ 



(1) 



where, e = 0.053, k = 107896Pa/m, ki = 9162Pa/m, 
K = 113919Pa/m for the dense foam and e = 0.0905, 
k = 39570Pa/m, ki = 3164Pa/m and K = 54704Pa/m 
for the light foam. As is expected the stress-strain be- 
havior is linear for strains less than a critical threshold: 
e < £q ~ 0.054 and e > ~ 0.09. For strains above the 
threshold the stress exhibits a large plateau as the foam 
cells undergo "collapse", signalled by a large change in 
strain at a nearly constant stress. Upon increasing the 
strain further the foam undergoes progressive densifica- 
tion as the opposite faces of foam cells meet, inducing 
the observed nonlinear stress-strain response. 

We model the foam as a homogeneous material made 
of nonlinear clastic elements obeying the force-extension 
relation given by Eq.([T]), translationally invariant in y 
and z directions. The foam is held fixed at x = and a 
stress of the form 
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FIG. 2: Strain dependence of 7, dense foam (+) top curve 
and light foam (x) bottom curve respectively. Solid lines are 
theoretical fits to the data. 



is applied at x = L. 

As shown in the supplementary material, the equation 
describing the non-linear foams considered here is given 
by 

du _ f u x > -e : k u xx . , 

1 dt \ u x < -e : kiu xx + K 6 (-e - u x ) b u xx 

where u(x) is the local displacement relative to the un- 
strained configuration of the foam; 7 is a friction coef- 
ficient which, a priori, can depend on the strain. We 
neglect the inertial term as in our set up it is 13 or- 
ders of magnitudes smaller than the damping term due 
to the low density of the foam. Notice that the bound- 
ary conditions above are such that u(0) = 0, u(x) < 
for compression and e = —u(L)/L. To determine 7, we 
have measured tan(<5) for the two foams using a dynamic 
strain of 0.001 and simulated these experiments using eq 
(|3"|. The strain dependent damping coefficient 7 was fit- 
ted to the experimental data using 



7 = 7o(l + 7i e") 



(4) 



with 70 = 5.31981 10 8 Pa/s 71 = 60 and n = 6 for the 
dense foam and 70 = 4.87 10 7 Pa/s 71 = 22.62 and n = 4 
for the light foam, with very good agreement as shown 
in figure O 

To study the properties of a foam subjected to a large 
amplitude periodic strain of period T, it is worth noting 
that eq ([3|) can be rescaled as follow: 



du 



u y > -e 



kg Uyy 



(1 + 76 u )~r = S u y < -e : k{u 



dr 



F(t) = F + F A cos[vt} 



(2) 



+ K* 6 (-60 - Uy) 5 Uyy 

(5) 

where, defining p = T/jqL 2 , y = x/L and r = t/T 
we have fcg = kop, k\ = kip, K* = Kp and F* = 
F p. The scaling parameter p allows us to relate the 
length scale, time scale and damping parameter 70 of 
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FIG. 3: Displacement u(x) of a foam subjected to periodic 
excitation F = 0.128075 * (1 + cos(t/T)) during 1 period of 
oscillation, fcj « 1.382, k{ « 0.1173, « 1.459 and 71 = 60. 
a) p = 0.01, b) p = 0.05. The sequence for the profiles is 
indexed in figure b. 



the foam. We thus see immediately that increasing the 
period of the load, T, decreasing the sample length L or 
the foam friction 7 are equivalent in this rescaled system 
and correspond to a variation of p. For the dense foam 
we have k* w 1.382, fc* w 0.1173, K* « 1.459 and F* w 
1.28075 10- 5 F/Pa 

The amount of energy transferred, -E tr , to the foam by 
the excitation can be easily evaluated as 



Ei 



T „du L f T du 

F—dt = - / F—dr. 
o dt p J dr 



(0) 



The energy dissipated by the foam, E^, on the other 
hand, is given by 



E 



dis 



JO 



l{ d 4Ydtd x = ± I" f T ( d ±) dr d y.{l) 



dt 



P Jo Jo 



dr 



(8) 



Then, the absorption quality of the foam will be given 
by the ratio Edi S / E tr which corresponds to the relative 
amount of absorbed energy. 

To analyse the properties of a foam under large ampli- 
tude periodic stress, we have considered the dense foam 
subjected to a sinusoidal pulse F(t) = 0.128*(l+cos(r)). 
The maximum stress F* = 0.256 correspond to a load of 
20kPa for the dense foam and as can be seen from fig [1] 
this correspond to a strain e = 0.6. We are thus probing 
the full range of the foam strain. 

On figure |3] we present snapshots of the foam displace- 
ment at regular interval during one period of oscillation, 
when subjected to the sinusoidal pulse described above. 
Fig [3]a corresponds to p = 0.01 while [3jb corresponds 
to p = 0.05. It is interesting to notice that the dis- 
placement profile is not sinusoidal as one might expect 
but instead progresses inside the foam as small plateaux. 
Small values of p correspond to large values of 7. One 
thus sees that in that case, the amplitude of excitations 
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FIG. 4: Strain (e = —u(L)/L) of a foam subjected to periodic 
excitation F = 0.128075 * (1 + cos(t/T)) during 1 period of 
oscillation. k* w 1.382, k* ~ 0.1173, K* « 1.459 and 71 = 60. 
p = 1.0, 0.3, 0.1, 0.07, 0.03, 0.01, 0.003, 0.001, 0.0001. 



are decreasing as one moves towards the origin. When 
p = 0.05, the foam deformation progresses as a pulse, 
exhibiting large regions where the strain, du/dx, is very 
small. At the base of the foam, one notices that the strain 
is very small for the first half of the cycle and then large 
and constant for the second half. 

On figure Q] we present the time evolution of the foam 
strain during one cycle. For large values of p, the strain 
oscillates between and 0.7 while for very small value 
it is nearly constant, just below 0.4. In the former case, 
the foam does not dissipate much and the load is trans- 
ferred nearly unaffected by the foam. In the later case, 
on the other hand, the foam is very dissipative and the 
foam adopts a constant strain effectively averaging over 
the load it is subjected to. It is also interesting to note 
that there is a range of values around p = 0.01 where 
the strain of the foam is nearly constant but with a very 
small value (« 0.25). This is the value of p such that 
the inside deformation of the foam are attenuated over 
a length scale comparable to the thickness of the foam. 
As a result, the foam is optimally absorbing. This is well 
illustrated on figure [5] where we show the relative absorp- 
tion of the foam as a function of p for various values of 
the nonlinear elasticity parameter K*. The maximum, 
which docs not vary much with K*, lies mostly in the 
range p = 0.01 to p = 0.001. When K* = 150 the foam 
is very stiff and the penetration much smaller, explaining 
the small absorption. Moreover, we can see from Fig- 
ure Ob that the phase difference between the excitation 
and the foam deformation is larger than 90 degrees in 
the region of maximum absorption. There is thus a par- 
tial phase opposition between the stress and the induced 
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FIG. 5: Dependence on p of the relative energy absorption 
(a) and Phase difference (b) of a foam subjected to periodic 
excitation F = 0.128075 * (1 + cos(t/T)). k* « 1.382, k* ~ 
0.1173 and 71 = 60. K* « 1.459 for (b) 



strain which results in the nearly constant low strain of 
the foam as described above. 

We have also considered a square pulse given by 
F{t) = 0.256(0 < r < t PW ) and F(t) = Q.12S{t PW < 
t < 1) and we have observed a similar behaviour of the 
foam. While of the oscillation profile inside the foam 
changed, the absorption properties of the foam was very 
similar to the one obtained for a sinusoidal excitation. 

In conclusion, we have derived an heuristic continuum 
model of dry foam deformation characterised by a non 
linear elastic response as well as a strain dependant fric- 
tion term. We have fitted the value of our model to the 
experimental data of two foams which both exhibited a 
6th order term for the inelastic term but a 6th and 4th 
power term for the strain dependence of the friction co- 
efficient. 

We have also studied the properties of the foam un- 
der large amplitude, periodic oscillation and have charac- 
terised a regime where the foam is optimally absorbing. 
The optimal range is frequency dependant and can be 
achieved by changing the thickness or the friction of the 
foam. This is potentially relevant for the design of pro- 
tective foam pads to protect fragile objects during trans- 
port. One should use a foam where the amplitude of the 
stress is such that it strains the foam in its non-linear 
regime. This will depend mostly on ko and K. If one can 
estimate the frequency of the excitations to which the 
object to protect will be subjected, for example 1Hz for 
a parcel being transported, then one can determine the 
optimal thickness of the protective pad after measuring 
7 for the selected foam. 
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